stable_log1p Function

public pure elemental function stable_log1p(x) result(rst)

Computes log(1 + x) accurately when x is near zero.

Arguments

Type IntentOptional Attributes Name
real(kind=real64), intent(in) :: x

Argument on [-1, infinity); NaN inputs propagate.

Return Value real(kind=real64)

Logarithm; x = -1 gives -infinity, x < -1 gives NaN.